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We investigate how emergent nematic order and nematic fluctuations affect several macroscopic 
properties of both the normal and superconducting states of the iron pnictides. Due to its magnetic 
origin, long-range nematic order enhances magnetic fluctuations, leaving distinctive signatures in 
the spin-lattice relaxation rate, the spin-spin correlation function, and the uniform magnetic sus- 
ceptibility. This enhancement of magnetic excitations is also manifested in the electronic spectral 
function, where a pseudogap can open at the hot spots of the Fermi surface. In the nematic phase, 
electrons are scattered by magnetic fluctuations that are anisotropic in momentum space, giving 
rise to a non-zero resistivity anisotropy whose sign changes between electron-doped and hole-doped 
compounds. We also show that due to the magneto-elastic coupling, nematic fluctuations soften 
the shear modulus in the normal state, but harden it in the superconducting state. The latter 
effect is an indirect consequence of the competition between magnetism and superconductivity, 
and also causes a suppression of the orthorhombic distortion below T c . We also demonstrate that 
ferro-orbital fluctuations enhance the nematic susceptibility, cooperatively promoting an electronic 
tetragonal symmetry-breaking. Finally, we argue that T c in the iron pnictides might be enhanced 
due to nematic fluctuations of magnetic origin. 



I. INTRODUCTION 

Understanding the normal state properties of the iron- 
based superconductors is a key step to decipher the na- 
ture of their high-temperature superconducting state (for 
reviews, see In fact, the superconducting (SC) tran- 
sition temperature T c is maximum near an instability to a 
spin-density wave state (SDW), suggesting that spin fluc- 
tuations may play an important role for the formation of 
the Cooper pairs Q. However, besides this SDW state, 
another ordered state is observed in the phase diagrams 
of the iron pnictides: the orthorhombic (Ort) phase. Un- 
like other superconductors that also display lattice in- 
stabilities - such as the cuprates - in the pnictides the 
Ort transition line always follows the SDW transition line 
across the phase diagram, even when the latter is bent- 
back inside the SC dome 0, 0| ■ Most importantly, the 
Ort transition at T s either precedes or is simultaneous 
to the SDW transition at T/v, implying that the lattice 
distortion is not a mere consequence of long-range mag- 
netic order. Since these two phases are closely related, 
and since magnetic fluctuations are the main candidate 
to explain the high-T c SC state, it is natural to investi- 
gate in more detail the origin of the Ort phase and its 
interplay with SC. 

On the experimental side, the development of detwin- 
ning techniques (see [f| for a review) allowed researchers 
to investigate the anisotropic properties of the Ort phase 
using several different experimental probes [54l3|. It 
became clear that the small lattice distortion could ex- 
plain neither the magnitude nor the doping-dependence 



of these anisotropics - in particular, the remarkable am- 
plitude of the resistivity anisotropy in electron-doped 
samples 0, [l4| and its sign-change upon hole doping [l5[ . 
Instead, this pool of observations suggest that the Ort 
phase is actually a manifestation of another electronic 
ordered phase that breaks the same tetragonal symme- 
try of the system. 

On the theoretical side, two different approaches have 
been proposed to explain the connection between the 
Ort and SDW phases, as well as the anisotropic prop- 
erties displayed below T s . In one approach, the onset of 
the SDW is regarded as a consequence of the Ort phase, 
which itself is driven by a spontaneous ferro-orbital or- 
der [l6l - [2l| . In this scenario, the anisotropic properties 
displayed in the Ort phase are a consequence of the un- 
equal occupations between the d xz and d yz orbitals of 
the Fe atom. Indeed, a spontaneous orbital polarization 
is obtained in some (but not all, see Ref. [ljl) strong- 
coupling models - in particular, in some versions of the 
Kugel-Khomskii model designed for the iron pnictides. 
However, in weak-coupling or moderate-coupling mod- 
els, orbital order is usually found only inside the SDW 
phase, but not as a spontaneous instability of the system 
preempting the SDW transition (2^ - [23 |. 

The other approach adopts the opposite point of view, 
i.e. that the SDW phase causes the Ort phase (25r - 
[35| . Since T s > Tjv, it is not long-range magnetic or- 
der, but magnetic fluctuations, that spontaneously break 
the tetragonal symmetry of the system. Interestingly, 
similar ideas were proposed in the early 1970s to ex- 
plain the splitting between the magnetic and the cubic- 
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Figure 1: Schematic representation of the nematic transi- 
tion in real space, (a) The transition from the disordered 
phase to the SDW phase breaks an 0(3) x Z2 symmetry. The 
0(3) symmetry refers to rotations in spin space while the Z2 
(Ising) symmetry refers to the two degenerate ground states 
of magnetic stripes with parallel spins along the y axis (order- 
ing vector Qi = (ir,0)) or along the x axis (ordering vector 
Q2 = (0, 7r)). (b) The 0(3) x Z2 symmetry can be broken in 
two steps. First, only the Z2 symmetry is broken: the system 
is still paramagnetic, since (Si) = (indicated by the gray 
double-arrow on top of the spins), but the spin correlations 
break the tetragonal symmetry, (Si ■ Si +X ) = — (Si • Si+ y ) 
(red and blue bonds, respectively). In the second step, the 
O (3) symmetry is broken and the system acquires long-range 
magnetic order. 



to-tetragonal transition in a family of rare-earth pnic- 
tides, such as DySb, CeSb, DyP, HoP, DyAs [H|. 

The qualitative idea is simple and can be understood 
using symmetry arguments, as shown in Fig. Q] In 
many antiferromagnets, the symmetry that is broken at 
the magnetic transition temperature is the 0(3) spin- 
rotational symmetry. To the 0(3) symmetry breaking 
corresponds also a translational symmetry breaking, due 
to the increase in the size of the crystalline unit cell in the 
magnetically ordered phase. In the iron pnictides, how- 
ever, the situation is different. The SDW ground state 
is actually doubly-degenerate, as it is characterized by 
magnetic stripes of parallel spins along either the y axis 
(ordering vector Qi = (jr, 0)) or the x axis (ordering vec- 



tor Q2 = (0,7r)). Therefore, to go to the ordered state, 
the system has to break not only the 0(3) spin-rotational 
symmetry, but it also has to choose between two degen- 
erate ground states, which corresponds to a Z2 (Ising- 
like) symmetry. Since Z2 is a discrete symmetry, the Z2 
symmetry-breaking is expected to be less affected by fluc- 
tuations than the continuous 0(3) symmetry-breaking, 
what opens up the possibility of the former happening 
before the latter. 

This is the idea behind the Ising-nematic state: an in- 
termediate phase preempting the SDW state, where the 
Z2 symmetry is broken but the 0(3) symmetry is not. In 
real space, the Z2 symmetry-breaking corresponds to a 
broken tetragonal symmetry, since the correlation func- 
tions (Si • Si +X ) and (Sj • Sj+y) acquire opposite signs 
(see Fig. DIb)). This is the origin of the term ne- 
matic: in liquid crystals, a nematic phase is character- 
ized by a broken rotational symmetry and an unbroken 
translational symmetry. Although the translational and 
rotational symmetries are always broken in crystalline 
solids, the analogy remains valid: in the electronic ne- 
matic phase, the point-group symmetry is reduced from 
04 (tetragonal) to 02 (orthorhombic), corresponding to 
an additional lowering of the rotational symmetry. From 
a purely symmetry point of view, the nematic state is 
therefore equivalent to the orthorhombic phase, which 
is the result of the inevitably induced distortion of the 
crystalline lattice (see below). The term "nematic" is 
however used to emphasize the fact that the phase tran- 
sition is of purely electronic origin and would still take 
place in a perfectly rigid lattice (for nematic transitions 
in other s yste ms, see Ref. (13|)- The recent analysis of 
Chu et aL[38l| demonstrated that it is indeed possible to 
experimentally determine the driving force behind the 
transition in the iron pnictides and distinguish between 
electronic driven transitions and ordinary tetragonal-to- 
orthogonal structural transitions. 

The mechanism behind the spontaneous breaking of 
the additional Ising symmetry is reminiscent of the 
order-out-of-disorder mechanism put forward by Chan- 
dra, Coleman, and Larkin in the context of localized 
spin models j39|. Not surprisingly, the first model cal- 
culations that obtained an spontaneous nematic phase 
in the pnictides were based on a strong-coupling ap- 
proach, the so-called J1-J2 model [2^, [26|, [28[. More 
recently, it was shown that an itinerant description of 
the system also accounts for a preemptive nematic phase 
[3o| . Notwithstanding important differences between the 
strong-coupling and weak-coupling approaches - in par- 
ticular on the character of the nematic and ma gne tic 
transitions as function of doping and pressure [35| - 
they share similar physics: magnetic fluctuations spon- 
taneously break the tetragonal symmetry already in the 
paramagnetic phase. 

In this paper, we will focus not on the origin of the 
nematic degrees of freedom, but rather on its manifesta- 
tions in several properties of the iron pnictides, such as 
the spin-spin correlation function, the spin-lattice relax- 
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Figure 2: Schematic phase diagram of the iron pnictides, em- 
phasizing the role played by nematic order and nematic fluc- 
tuations. SDW denotes the spin-density wave state, SC the 
superconducting phase, PM the paramagnetic phase, and Tet 
the tetragonal phase. Tetragonal symmetry is only broken 
below the nematic/orthorhombic transition line, but nematic 
fluctuations remain at higher temperatures, as evidence d by 
shear modulus [UEOl, uniform susceptibility anisotropy 
and resistivity anisotropy 0, 0] measurements. 



manifested in the suppression of the orthorhombic distor- 
tion below T c and in the hardening of the shear modulus 
inside the SC dome. Despite this competition, we will 
give a qualitative argument of why nematic fluctuations 
can enhance the value of T c without affecting the symme- 
try of the SC state. The schematic phase diagram shown 
in Fig. [2] summarizes the main points of this work: not 
only nematic order, but also nematic fluctuations, play 
an important role in the physics of the iron pnictides. 

The paper is organized as follows: in Section (TT] we re- 
view the model that relates the nematic order parameter 
and the nematic susceptibility to the magnetic spectrum 
of the system. In Section Hill we use this model to calcu- 
late several system properties that can be probed exper- 
imentally: the magnetic properties in the PM-Ort phase 
(spin-lattice relaxation rate, spin-spin correlation func- 
tion, and uniform susceptibility); the electronic spectral 
function in the PM-Ort phase; the orthorhombic distor- 
tion in the PM-Ort phase and the shear modulus in the 
PM-Tet phase; and the resistivity anisotropy above the 
SDW transition temperature. We address the interplay 
between nematic and ferro-orbital degrees of freedom in 
Section Hvl while the interplay with superconductivity is 
investigated in Section [Vj Section I VII is devoted to the 
conclusions and closing remarks. 



ation rate, the resistivity anisotropy, the uniform mag- 
netic susceptibility, and the elastic modulus. Since all 
these quantities are experimentally accessible, they pro- 
vide important cornerstones to test the predictions of the 
nematic model. In particular, we will investigate not only 
the role played by nematic order but also by nematic 
fluctuations. The former is trivially manifested as a fi- 
nite orthorhombic distortion of the lattice, but, due to 
its magnetic origin, it also leaves distinctive signatures 
in several magnetic properties. For instance, the onset 
of nematic order enhances magnetic fluctuations, caus- 
ing a kink in the spin-lattice relaxation rate, and splits 
the degenerate spectrum of magnetic excitations around 
Qi = (7r,0) and Q2 = (0,7r). Furthermore, we will show 
that the anisotropic uniform susceptibility and the resis- 
tivity anisotropy are proportional to the nematic order 
parameter near T s , although the sign of the latter de- 
pends on the geometry of the Fermi surface. 

Besides nematic order, low-energy nematic fluctua- 
tions also emerge, strongly affecting the system proper- 
ties. These fluctuations renormalize the shear modulus 
already at high temperatures, providing a way of prob- 
ing the nematic susceptibility directly from elastic mea- 
surements. We will also discuss the interplay between 
nematicity and the other electronic degrees of freedom of 
the system. In particular, we will show how ferro-orbital 
fluctuations can enhance the tendency to a nematic insta- 
bility, which in turn can promote orbital order. Finally, 
the interaction between nematic and superconducting de- 
grees of freedom will be addressed. As a consequence of 
its magnetic origin and of the competition between SDW 
and SC, nematicity indirectly competes with SC. This is 



II. MAGNETIC MODEL FOR THE NEMATIC 
PHASE 

We start with an effective action describing the mag- 
netic degrees of freedom (35L l40j . Since there are two 
magnetic ground states, we introduce two order param- 
eters Mi and M2 corresponding respectively to stripes 
with parallel spins along the y axis (i.e. ordering vec- 
tor Qi = (-7T, 0)) and stripes with parallel spins along 
the x axis (i.e. ordering vector Q2 = (0,7r)). Thus, 
the spatial magnetic configuration is given by S (r) = 
M l6 »Qi-r + M 2 e lC * 2 ' r . The free energy describing these 
degrees of freedom is given by: 

-Fmag = / Xo 1 (q) (Mi,, • Ml_, + M 2 , 9 • M 2 ,_,) 

Jq 

+ 1 jf (Mi + Mlf - I jf {M\ - Mlf (1) 

where J q = Tj^nl ikh and 9 = (q> w n) denotes the 
momentum q and the Matsubara frequency <j„ = 2nnT. 
Here x^ 1 (q) = tq + q' 2 + 7 _1 \uj n \ is the bare dynamic 
spin susceptibility and tq oc {T — Tjv,o); with Ijv^o the 
mean-field SDW transition temperature and 7 the Lan- 
dau damping. Notice that M^ g is a function of both mo- 
mentum and Matsubara frequency, as it has both static 
and dynamic fluctuations. 

The coupling constants u and g determine the nature 
of the magnetic ground state. Consider, for instance, 
the uniform limit and minimize the free energy F with 
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respect to Mj. Then, if g > and u > 0, we obtain 
two solutions for T < Xjv,o corresponding to M 2 = and 
Mi = yj—ro/u or M\ = and M2 = \/—ro/u. These are 
nothing but the two degenerate magnetic stripe configu- 
rations we discussed earlier. Therefore, hereafter we con- 
sider u > and g > 0. From symmetry considerations, 
the action (JlJ can also contain the term (Mi • M2) , 
which, for g > 0, does not affect the selection of the 
ground state as long as its coefficient is not a large neg- 
ative number compared to g. 

Although here we introduced the free energy ([1]) in a 
phenomenological way, it can be derived microscopically 
from both the itinerant [35[ and localized-spin approaches 
[25| . In both cases, the term (Mi • M2) 2 is absent. Even 
though both approaches find g > 0, which selects the 
magnetic ground states observed in the iron pnictides, 
the magnitude of this coupling constant is different in 
the band model and in the J\ — J2 model. In particular, 
while in the latter g is small and induced by thermal or 
quantum fluctuations, in the former it can be large, as 
it is given by products of the non-interacting electronics 
Green's functions. It has been suggested that a relatively 
large g is necessary in order to fit the spin- wave spectrum 
inside the SDW phase [3l| . 

In terms of the magnetic degrees of freedom, the ne- 
matic order parameter can be expressed as ip oc M^—M^. 
Indeed, (ip) 7^ implies that the fluctuations around 
one of the ordering vectors (Mf\ are on average dif- 
ferent than the fluctuations around the other ordering 
vector (Mf). Since Qi = (n, 0) and Q 2 = (0, n) are 
related by a 90° rotation in the square-lattice unit cell, 
(Mi ) 7^ (M'2) implies that the x and y directions are in- 
equivalent and the tetragonal symmetry is broken. This 
tetragonal symmetry-breaking is enforced by magnetic 
fluctuations, and not long-range magnetic order, since 
(Mi ) ^ (M|) does not require or implies (M;) = 0. 

Once long-range nematic order sets in, it has a feed- 
back effect on the magnetic fluctuations, assisting the 
transition to the SDW state. Physically, the nematic or- 
der parameter "transfers" magnetic spectral weight from 
one ordering vector to the other, enhancing magnetic 
fluctuations in one channel but suppressing them on the 
other (see Section UnC ). This process increases the SDW 
transition temperature, as it becomes larger than what 
it would be in the absence of long-range nematic order. 

The expressions relating the nematic order parame- 
ter ip and the magnetic susceptibility can be obtained 
by introducing two auxiliary Hubbard-Stratonovich fields 
in the free energy ffl, one for each quartic term (for 
more details, see J35j]). The first one accounts for 
the Gaussian fluctuations of the magnetic order pa- 
rameter, (M 2 + M|), whereas the second one accounts 
for the possibility of a finite nematic order parameter, 
(Mi — Mf). In the saddle-point (1/^V) approximation 
[4l(, one obtains a set of non-linear coupled equations 
for these two auxiliary fields: 



■•P 



77 / ixi (?) - X2 (<?)] 



r = r + ~ [xi (q) + Xi (<?)] 



(2) 



where r oc £ -2 , with £ denoting the magnetic correlation 
length renormalized by Gaussian fluctuations. Xi (?) i s 
the renormalized magnetic susceptibility at the ordering 
vector Q.;: 



xi (q) 

X2 (q) 



1 



(r -tp) + q 2 + 7- 1 |w„| 
1 

(r + ^) + g 2 + 7- 1 |w„| 



(3) 



Expanding equations ([2]) for small tp, we obtain a fa- 
miliar </> 4 -type equation of state for tp: 



(4) 



1 - a / X 2 (q) )<P + bp 3 + O {p 5 ) = 



where % (q) has the same form as the bare susceptibility 
Xo (q) but with ro replaced by r, i.e. x~ l (?) = ?' + 
g 2 + 7~ 1 |w„|. The sign of the cubic coefficient b depends 
on the ratio u/g and on the dimensionality d (35[, and 
determines whether the nematic transition is first-order 
or second-order. 

For b > 0, the sign of the linear coefficient determines 
when a finite ip 7^ becomes a solution of the equation 
of state (TJJ. At high enough temperatures, the magnetic 
fluctuations are weak and J x 2 (?) x d 1S small, im- 
plying that the linear coefficient is positive. Near a mag- 
netic instability, however, f x 2 (?) increases significantly, 
since it diverges at the magnetic transition temperature 
for d+z < 4 (z is the dynamic critical exponent, relevant 
for quantum phase transitions). Therefore, proximity to 
an SDW transition makes the linear coefficient change 
sign, inducing the nematic phase transition. Notice that, 
for b < 0, the first-order transition may take place while 
the linear coefficient is still positive. 

Another way of expressing the same results is via the 
nematic susceptibility, defined as the correlation function 

(q) = ((A/ x 2 - Mf) 9 (Ml - A/, 2 ) _V After intro- 
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free energy ([T]) and then taking h — > 0, we obtain: 



ducing a symmetry-breaking term h (Ml — M% ) in the 
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In the saddle-point approximation, it follows that [21 
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f q X 2 (q) 

1 - 9 L x 2 (q) 
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Figure 3: Diagrammatic representation of the nematic sus- 
ceptibility in terms of the electronic operators Ck, /i,k = 
Ck+Qi +c k +Q 2 and f 2 ,k = Ck+ch -Ck+q 2 , where c k is the oper- 
ator associated with states around the center of the Brillouin 
zone and Qi are the magnetic ordering vectors. Dashed lines 
are the external nematic legs, solid lines are the electronic 
propagators, and the wavy lines are the magnetic propaga- 
tors associated with Li = Mi + M2 and L2 = Mi — M2. 
The nematic susceptibility (on the right) corresponds to an 
Aslamazov-Larkin type diagram. 



As expected, the condition for the onset of long-range 
magnetic order, x ncm — > 00, agrees with the condition for 
a finite tp ^ in the equation of state (U]). 

It is interesting to obtain a more microscopic under- 
standing of the nematic susceptibility in terms of the 
electrons. For this, we consider the states around the 
center of the square-lattice Brillouin zone (annihilation 
operator Ck), where the hole pockets are present, and the 
states Ck+Qi in the electron pockets centered at the the 
magnetic ordering vectors Qi = (tt, 0) and Q2 = (0, 7r). 
The basic nematic vertex is given by ip (Mj 2 — Mf ) , 
and the collective spin modes are expressed in terms 
of the electronic operators as Mj = ^ k c^ a tr a pCk+Q ( f), 
with Pauli matrices cr a p. It is convenient to intro- 
duce the two-sublattice Neel vectors Li and L 2 , de- 
fined as Li = Mi + M 2 and L 2 = Mi - M 2 . Then, 
it follows that Lj = c ka <r a| g/i,k/3, where we intro- 
duced the fermionic operators /i k = Ck+Qj + Ck+Q 2 and 
/ 2 ,k = Ck+Qj — Ck+Q 2 . Within this notation, the nematic 
vertex becomes ip (Li • L 2 ), and the nematic susceptibil- 
ity is given by an Aslamazov-Larkin type diagram, as 
shown schematically in Fig. [3] 

Notice that Xncm is a four-spin correlation function. 
In a regular antiferromagnet, where only one order pa- 
rameter is present, its critical behavior can be fully ex- 
pressed in terms of the critical behavior of the two-spin 
correlation function - the magnetic susceptibility 
In that case, the four-spin correlation function does not 
contain any information that is not already embedded in 
the two-spin function. In the present case, however, this 
is not true: due to the nematic degrees of freedom, the 
four-spin correlation function acquires its own critical be- 
havior, revealing the emergent character of the nematic 
phase. 

In our subsequent analysis we primarily focus on the 
anomalies that take place at the nematic transition. For 
a realistic description of the pnictides the behavior at 
the nearby magnetic transition is of course closely con- 
nected to that at the Neel temperature (see Fig. HJ. 
However, in order to stress the behavior of the param- 
agnetic nematic phase we will frequently consider either 



two-dimensional systems, where the magnetic ordering 
temperature is suppressed to zero, or systems where the 
splitting between the two phase transitions is not too 
small. 



III. MANIFESTATIONS OF NEMATIC ORDER 
AND NEMATIC FLUCTUATIONS IN THE 
NORMAL STATE 

One immediate consequence of the coupled non-linear 
equations @ is that the nematic and magnetic transi- 
tions tend to follow each other. For instance, the di- 
vergence of the static nematic susceptibility in Eq. © 
depends on how close the system is to a magnetic insta- 
bility, where J x 2 (?) diverges. At the same time, once 
long-range nematic order sets in, the static magnetic sus- 
ceptibility Xi (<7 = 0) is enhanced, as shown in Eq. ([3]), 
bringing the SDW transition temperature up. Therefore, 
both the nematic and the magnetic transitions are tied 
together, in qualitative agreement with the experimental 
phase diagrams of the iron pnictides. This is perhaps the 
most striking feature in favor of the nematic model. 

The behavior of the solutions of Eq. ([2]) has been stud- 
ied in detail elsewhere, and the characters of the magnetic 
and nematic transitions have been investigated for a wide 
range of parameters [35j . In the remainder of the paper, 
we will focus instead on the manifestations of long-range 
nematic order, as well as nematic fluctuations, in several 
properties of the system. 

A. Orthorhombic distortion: x-ray diffraction 

In the nematic model, the spontaneous breaking of the 
tetragonal symmetry takes place in the magnetic sector, 
since magnetic fluctuations become stronger around one 
of the two ordering vectors Qi = (tt, 0) and Q 2 = (0, tt), 
which are related by a 90° rotation. Following Landau's 
theory of phase transitions, once the tetragonal symme- 
try is broken in the magnetic sector, it will be broken in 
all sectors at the same temperature. Indeed, if one con- 
structs another (scalar) order parameter rj that breaks 
tetragonal symmetry, it must couple bi-linearly to the 
nematic order parameter ip in the free-energy expansion, 
implying that rj cx p in the free-energy minimum. 

This is the case for the orthorhombic order parameter 
e s = (a — b) I (a + 6), where a and b are the lattice con- 
stants. Considering a harmonic lattice, the elastic free 
energy is given by: 

Fcl = / T £ " " Acl / £s " M ^ (7) 

where C s is the shear modulus and A e i is the magneto- 
elastic coupling. If one works in the coordinate system 
of the square lattice with one Fe per unit cell, then C s = 
C\\ — C\2\ instead, if one works in the coordinate system 
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Figure 4: Nematic order parameter (p (in units of its value 
ifN at the magnetic transition temperature Tjv) as function 
of the reduced temperature T 3 —T (arbitrary units). In panel 
(a), the split magnetic transition is first-order, while in (b) 
and (c) it is second-order and progressively farther from the 
magnetic tricritical point, ip is obtained by solving Eqs. ([2]) 
for a three-dimensional anisotropic magnetic dispersion. The 
only parameter changed in each plot is the ratio u/g. 



ter {Ml) - (Ml) in Eq. Q, see Ref. [H|. Then > onc 
has to add to the elastic free energy the term ae s : 

F ci = f ^ e 2 s - A e i f e s {Ml - M|) - f ae s (8) 

J X * J x J X 

For a harmonic lattice, the elastic contribution to the 
partition function Z — J dMi de s e~ Fn,a!S ~ Fo1 can be inte- 
grated out analytically, yielding the renormalized mag- 
netic free energy: 




Comparing to the original action, Eq. the nematic 
coupling g is enhanced (2tJ (which holds also for twin 
samples) and a is converted into an external field for the 
nematic order parameter, whose amplitude depends on 
A e i/C s . The impact of this external field to the nematic 
and magnetic transition temperatures has been widely 
discussed both theoretically Q |45[ and experimentally 

SHI, S3. 

B. NMR \jT\T spin-lattice relaxation rate 



with two Fe per unit cell, then C s = Cqq. In either case, 
minimization of Eq. §7$ leads to e s cx (Ml) — (M|), 
i.e. measuring the orthorhombic distortion is equivalent 
to measuring the nematic order parameter. Of course, 
the mere detection of an orthorhombic distortion is by 
no means an evidence for the electronic character of the 
structural transition. Yet, one can ask whether there are 
any signatures of the magnetic character of the structural 
transition in the behavior of s s . 

Interestingly, the solution of the set of self-consistent 
equations @ shows that, when the magnetic and ne- 
matic transitions are split, e s is in general not affected 
by the magnetic transition, except near a magnetic tri- 
critical point of the phase diagram [35j . In this re- 
gion of parameter-space, e s has a kink at the magnetic 
transition temperature - more specifically, \de s / dT\ T - < 
\de s /dT\ T + , as shown in Fig. 0] Away from the tri- 
critical point, however, e s becomes a smoother function 
around T/y. Indeed, x-ray diffraction measurements in 
Ba (Fei_ x Co x ) 2 As2 find a kink in the orthorhombic dis- 
tortion only in the proximity of x ~ 0.02 [42j, where a 
magnetic tricritical point has been experimentally estab- 
lished [H . 

Notice also that the magneto-elastic coupling has im- 
portant implications for detwinned samples. Detwin- 
ning can be achieved by applying a strain e s , which 
acts as an external field to the nematic order parame- 



In the nematic model, magnetic fluctuations are re- 
sponsible for the structural transition. Thus, it is natural 
to expect that the magnetic properties of the system will 
undergo noticeable changes at the structural (nematic) 
transition temperature T s . We first study the 1/T\T 
spin-lattice relaxation rate, which can be directly probed 
by NMR: 

q 

Here 7 g is the gyromagnetic ratio and A (q) is the 
structure factor of the hyperfine interaction. In our 
system, Xmag (q, ^) is strongly peaked at the order- 
ing vectors Qi = (7r,0) and Q2 = (0, 7r), with 
Xmag (q + Qi, w) = Xi (q, w n -> -icu + 0+) given by Eqs. 
([5]). Assuming A (q) constant, and considering a two- 
dimensional system for simplicity, we can evaluate the 
momentum sum around each ordering vector and obtain: 

where n > is a constant. In the absence of nematic 
order, (p = 0, one recovers the usual result 1/T\T cx 
1/ (T — Tjy). However, if the paramagnetic phase under- 
goes a nematic transition, the spin-lattice relaxation rate 
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Figure 5: Inverse of the spin-lattice relaxation rate T\T 
as function of the reduced temperature T — T a for a two- 
dimensional system. The change in the slope at T s is evident, 
and should be present whenever the magnetic and structural 
transitions are split. 



diverges when r = \ip\, which is the condition for the 
magnetic transition to take place, see Eq. (|5|). 

Most interestingly, when the nematic/structural and 
magnetic transitions are split, 1/T\T displays a kink at 
the nematic transition, since: 



1 



1 



oc 



(12) 



In Fig. [5l we illustrate this behavior by presenting 
the temperature dependence of T±T as calculated from 
the solution of Eq. ([2]) at d = 2. The change of slope 
at T s is evident. Experimentally, this feature should be 
more pronounced and easier to observe in compounds 
with well separated magnetic and structural transitions. 
Indeed, in NaFeAs, which has T N w 40 K and T s w 60 
K, NMR measurements (48l |49| have found a very clear 
change in the slope of 1/T\T at T s , in agreement with 
the predictions of the nematic model. 



C. Spin-spin correlation function: neutron 
scattering 

A more direct probe of the impact of long-range ne- 
matic order on the magnetic spectrum can be given 
by inelastic neutron scattering, as it directly mea- 
sures the spin-spin correlation function 5 (q, w) = 
Im [xmag (q 5 ui)] /uj at the two ordering vectors Qi = 
(tt, 0) and Q2 = (0,7r). In the paramagnetic phase, we 
have the usual expressions for overdamped spin excita- 
tions: 



S (q+Qi,w) 



So 



[(r -<p) + (l + ri) q 2 x + (1 - rj) ql] " + 7~ 2 ^ 2 



(13) 



| 4 
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Figure 6: (a) Spin-spin correlation functions Si (q, 0) = 
S (q x + Qi, w) (blue) and S2 (0, q) = S (q y + Q2, ui) (orange) 
as function of momentum transfer q for fixed energy cj/j — 0.1 
and several temperatures. The upper pairs of curves (coinci- 
dent) are for a temperature immediately above the structural 
transition T s while the lower pairs of curves refer to temper- 
atures progressively smaller than T 3 . Here the calculations 
were done at d = 2. In (b) and (c), as function of the reduced 
temperature T — T s , we show the q = peaks and the half- 
width at half-maximum HWHM of each correlation function 
51 (<?,0) (blue) and S2 (0,q) (orange). 



and: 



S (q+ Qa,w) 



So 



[(r + ip) + (l-ri)qZ + (l + ri)q*\ + 7~ 2 w 2 



(14) 



with r cx £ -2 and ip solutions of the self-consistent equa- 
tions @ and — 1 < 77 < 1 . These expressions were derived 
from Eq. §3§ with the isotropic magnetic dispersion q 2 re- 
placed by the anisotropic dispersion (1 ± 77) g 2 + (l T v) Qy 
. to make comparison with experiments more realistic 

The nematic order parameter can be obtained by 
comparing the correlation functions S (q x + Qi, u>) and 
S (q y + Q2,w) below T s but still above Tjv (i.e. in the 
paramagnetic phase). For a fixed low energy w/7 C r, 
they have peaks at q = proportional to (r ± ip)" 1 and 

half-widths proportional to (r ± ip)^ 1 ^ 2 . In Fig. [6l we 
plot the temperature evolution of both S (q x + Qx, ui) 
and S (q y + Q2, w), as well as the temperature depen- 



8 



dence of their peaks and widths. We used the solutions 
of the self-consistent equations © for d = 2, but the 
features at T s are of course much more general. 

Physically, the onset of long-range nematic order in the 
paramagnetic phase enhances the magnetic fluctuations 
around one of the two ordering vectors (in this example, 
Qi) while suppressing the fluctuations around the other 
ordering vector (Qa)- As a result, the constant-energy 
cuts of <S(q + Qi,w) become narrower and stronger, 
while the constant-energy cuts of 6> (q+ Q2,w) become 
broader and weaker. 

A neutron scattering experiment that detects this be- 
havior will provide perhaps the most direct evidence for 
the magnetic origin of the tetragonal symmetry-breaking 
in the iron pnictides. Of course, one needs a single crystal 
that displays split magnetic and structural transitions, 
such as the 1111 and the 111 compounds, or even some 
underdoped 122 compounds. The main experimental dif- 
ficulty is that the crystal needs to be detwinned - i.e. it 
must have only one structural domain. With the im- 
provement of detwinning techniques, it is reasonable to 
expect that such an experimental set-up will be available 
soon. 



D. Pseudogap in the spectral function: ARPES 

In the previous two subsections, we discussed the ef- 
fects of nematic order on the spectrum of magnetic ex- 
citations, and how these effects can be directly probed 
via magnetic measurements. Another alternative is to 
probe them indirectly via the electronic spectral func- 
tion, A(q, uS) = — 21m [G (q, uj)}, where G (q, w) is the 
electronic Green's function. 

Below Tjv, the electronic band structure is recon- 
structed by long-range magnetic order, and A (q, w) is 
modified due to the folding of shadow bands and the 
opening of a SDW gap. Above T^r, in the paramag- 
netic phase, the electronic spectral function is affected by 
magnetic fluctuations. In particular, several works have 
shown that magnetic precursors are able to reduce the 
spectral weight of the hot spots at the Fermi level (52l |. 
The position q^ s of these hot spots are given by the re- 
lation E (qhs) = E (q/j S + Qi), where E (q) denotes the 
band dispersion. To find the hot spots geometrically, one 
simply makes a copy of the Fermi surface and displaces 
it by the magnetic ordering vector Q^: the points that 
overlap with the original Fermi surface correspond to the 
hot spots. 

Deep inside the magnetically ordered phase, the spec- 
tral function at the hot spots vanishes at the Fermi 
level, i.e. A(qhs,oj = 0) = 0. At very high tempera- 
tures, A{(\h s ,uj) is maximum at u> = 0. At intermedi- 
ate temperatures in the paramagnetic phase, although 
A(c\j is ,uj = 0) 7^ always, strong magnetic fluctuations 
are able to suppress the zero-energy spectral weight at 
(\h s . Above a certain threshold of the magnetic correla- 
tion length £, the maximum at A(c\j is ,uj — 0) becomes 



A(q hs ,co) A(q hs ,co) 




T > T s T < T s 

Figure 7: Schematic representation of the opening of a pseu- 
dogap in the spectral function A (qh s ,a>) at the hot spot q/, s 
due to long-range nematic order. At high temperatures, above 
the nematic transition temperature T s , the magnetic correla- 
tion length is below a threshold £* and the peak of A (qhs,oj) 
is at uj — 0. The onset of nematic order at T 3 sharply enhances 
the magnetic fluctuations, as it was shown in Fig. [SJ and £ 
can overcome £*. In this case, the maximum of A(qhs,w) 
shifts to u « Asdw, which is the SDW gap at T — 0, and 
ui = becomes a minimum. 

a minimum as a result of spectral weight being trans- 
ferred to energies comparable to the zero-temperature 
SDW gap Agow j i- e - the maximum is displaced to 
A(qf ls ,u) = Asdw), see Fig. [7] Consequently, the elec- 
tronic spectrum has a pseudogap due to the presence of 
strong magnetic precursors. 

We can now understand how nematic order can affect 
the electronic spectrum. By sharply enhancing the mag- 
netic fluctuations, as it was shown before in Fig. [SJ the 
onset of long-range nematic order may trigger the open- 
ing of a pseudogap at the hot spots of the Fermi surface 
[35| . In fact, several experiments on the iron pnictides 
have reported possible manifestations of a pseudogap be- 
havior |53l . HH , although more work is still necessary to 
pinpoint its exact temperature scale and its relation to 
the magnetic and superconducting energy scales. The 
mechanism proposed here for the pseudogap could be di- 
rectly probed by ARPES or STM measurements, which 
are sensitive to the electronic spectral function. Probably 
the most promising material is NaFeAs, since T s and Tjv 
are already below 100 K and separated by approximately 
20 K. 



E. Uniform susceptibility: torque magnetometry 

Besides affecting the magnetic fluctuations at Qi = 
(■7T, 0) and Q2 = (0, it), long-range nematic order also 
triggers an anisotropy in the uniform magnetic suscepti- 
bility Xmag (q = 0). To calculate it, we include the contri- 
bution of a non-zero magnetic field H to the free energy 
((T|) and compare the effects of a field applied along the x 
direction, H x , and a field applied along the y direction, 
H y . We can either use symmetry considerations or per- 
form an explicit calculation by adding the Zeeman term 
to the electronic dispersions. Neglecting the cross terms 
H x H yi we obtain the general form: 
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F = F mag [Mi] + a x (H 2 X + Hi) (Ml + M, 



(15) 



+a 2 [Hi {Ml, + Ml,) + Hi (Ml y + M 2 2 , y ) ] 

with coupling constants a\ and 012. The first term is 
isotropic and describes the suppression of SDW order in 
the presence of an external field. The second term gives 
an anisotropic response. Since Xmag (q = 0) oc cPF/dHj, 
we obtain for the anisotropic uniform magnetic suscepti- 
bility 



{Ml x + Ml 



o)-xZ g (q = o) « 

' *'* )-(Ml v + Ml y ) 



(16) 



To understand the anisotropy in the susceptibility, we 
cannot use the 1/N (saddle-point) approach discussed in 
Section HH since it does not distinguish between the com- 
ponents of the magnetic order parameter. Yet, we can 
gain qualitative understanding by noting that, in the iron 
pnictides, the magnetic order parameter points parallel to 
the modulation axis in the SDW phase, i.e. Mi || x while 
M2 |j y. Therefore, at least within a phenomenological 
approach, there must be a spin-anisotropy term in the 
free energy expansion (JTJ) favoring these spin directions. 
Consequently, the main contribution to the nematic order 
parameter ip = Ml - Ml comes from (Ml x ) - (M% y ), 
implying that, to leading order in the spin-anisotropy: 



Imag 



(q = 0) - xSL (q = 0) oc ^ 



(17) 



Interestingly, recent torque magnetometry experi- 
ments on BaFe2As2 and BaFe2 (Asi_ x P x ) 2 measured the 
anisotropy in the uniform magnetic susceptibility, finding 
a non-zero value below the structural transition temper- 
ature fl3| . 



F. Resistivity anisotropy 



temperature T s . Of course, below Tjv the anisotropic 
reconstruction of the Fermi surface by long-range mag- 
netic order will give another contribution for Ap |22,[24|. 
but here we will focus only on the paramagnetic nematic 
phase. 

The electrons that are most efficiently scattered by 
spin fluctuations are the ones near the hot spots of the 
Fermi surface, E (q^s) = E (q/j S + Qi) ; since they re- 
main near the Fermi level even after acquiring the addi- 
tional momentum Q; transferred in the scattering pro- 
cess. The electrons on the other parts of the Fermi sur- 
face (called cold regions) do not significantly contribute 
to the spin-fluctuation scattering. For a very clean sys- 
tem, it was shown in Ref. that these cold regions 
of the Fermi surface dominate the transport properties, 
short-circuiting the contribution from the hot spots. On 
the other hand, Ref. [57} showed that in dirty sys- 
tems, the contribution of the spin-fluctuation scattering 
is the leading-order correction to the residual resistivity 
po. Applying these results to our system, since the re- 
sistivity anisotropy comes from the contribution of the 
hot spots, we expect Ap to be larger in samples that are 
dirtier. In fact, this has been recently confirmed by mea- 
surements comparing the resistivity anisotropy of normal 
and annealed samples (lp| . 

To proceed, we consider the dirty limit where po is 
large compared to the contribution of the spin- fluctuation 
scattering. We also consider a three-band model, with a 
central circular hole pocket T and two elliptical electron 
pockets X and Y displaced from the center by the or- 
dering vectors Qi = (n,0) and Q 2 = (0, 71"), respectively. 
Due to the tetragonal symmetry of the system, X trans- 
forms to Y under a 90° rotation. Within this model, the 
resistivity along the j axis is given by (see Ref. (55| for 
more details): 



Pjj - Po 
Po 



E 

a=X,Y 



; r, q 



v U) , 

a,q' 



,(«) 



,(<0 



+ t 



(18) 



The experimental tool that has been mostly used to 
probe the nematic phase is the resistivity anisotropy 
Ap = pb — p a d, 0, EL EH]. On the theory side, however, 
unlike the other properties discussed above, Ap refers 
to a non-equilibrium situation, what brings additional 
difficulties to the calculation. Here we follow the same 
framework of Ref. [55[ and employ a semi-classical Boltz- 
mann equation approach to construct a transport theory 
for the nematic phase. In this subsection, b and a are 
parallel to the y and x axis, respectively. 

The main property of the nematic phase is that mag- 
netic fluctuations are different around the two ordering 
vectors Qi = (n,0) and Q 2 = (0,7r). The scattering of 
electrons by these anisotropic spin fluctuations will give 
rise to an anisotropic scattering rate, resulting in a non- 
zero resistivity anisotropy Ap above the magnetic transi- 
tion temperature Tjv but below the structural transition 



where t ~ T/7 is a dimensionless temperature, vj)^ is 
the j-component of the Fermi velocity of band a, and 
c > is a constant proportional to the ratio between the 
amplitudes of impurity scattering and spin-fluctuation 
scattering. Here, q and q' are momenta restricted to 
the Fermi surface. The functions ^q^q' can be written 

as w q a j q , = (r ± tp) + |q r - q^ - Q a | 2 , where the upper 
(lower) sign refers to band Y (X). 

Performing a 90° rotation, and using the symmetry 
properties relating Vx to vy, we can expand Eq. ([T| 
for small (p near the structural transition, obtaining: 



Pyy Pxx 

Po 



tJ (y) 



c< ^E 

q.q' 



x/(qr-qx - Qi) 



(' 



(x) 



» 



r,q V X,q> 



(19) 
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Figure 8: Position of the hot spots in hole-doped and electron- 
doped compounds. The elliptical electron pocket (red) was 
displaced by the ordering vector Qi = (n, 0), crossing the cir- 
cular hole pocket (blue) at the hot spots. The changes in the 
position of the hot spots are a consequence of the changes in 
the areas of the electron and hole pockets due to the variation 
in the charge carrier concentration. We also present the Fermi 
velocity vp = — Vr,q hB at the hot spot in each case. Notice 
that Up > Vp in the hole-doped case, whereas > Up in 
the electron-doped case. According to Eq. (|19|) . one expects 
Ap < and Ap > 0, respectively. 



where the function / (k) is given by: 



/(k) 



2 (2k 2 + 2r + t) 



(k 2 



ty 



(20) 



Therefore, as expected, Ap ~ ip, but the sign of 
the proportionality constant can be positive or negative. 
Also, Ap may be small even when ip is large (or vice- 
versa) due to the interplay between impurity-scattering 
and spin-fluctuation scattering. 

The function / (fc) is in general dominated by the 
k = contribution - this can be rigorously shown at low 
temperatures and near an SDW quantum critical point, 
where r/t — > 0. Thus, the main contribution to the mo- 
mentum sum comes from q, q' near q/^q/is + Qi, and 
one can obtain the sign of the resistivity anisotropy with 
respect to <p by just analyzing the Fermi velocities of the 
hot spots. 

For instance, if the hot spots are close to the y axis, 
as it is the case for electron-doped samples (58|, then 
one expects the y-component of the Fermi velocity to 
be larger than the ^-component, i.e. Ap > 0. On the 
other hand, if the hot spots are close to the x axis, as 
it is the case for hole-doped samples, the opposite holds 
and one expects Ap < (see Fig. [5]). Hence, this general 
model predicts different signs of Ap (in the paramagnetic 
phase) for electron-doped and hole-doped systems. This 
feature has been recently observed experimentally [l5j . 
suggesting that the transport properties in the nematic 
phase are well described by this model of anisotropic spin- 
fluctuation scattering. 

We note that there are other contributions that 
may affect the resistivity anisotropy in the nematic- 
paramagnetic phase. For instance, orbital order (see 
Section IIV[) . which is also induced by nematic order, 
changes the Drude weight along the x and y directions. 
Whether this Drude weight redistribution can account 



for the experimentally observed resistivity anisotropy re- 
mains an unsettled issue, as some works find that the 
resistivity anisotropy coming from this mechanism alone 
gives the opp osite sign of the Ap measured experimen- 
tally [HI, [2C| Impurity effects also play an impor- 
tant role in transport, as discussed above and also in a 
different model proposed in Ref. [78[ based on orbital 
order induced by antiferro-orbital fluctuations. A recent 
Monte Carlo simulation in the spin-fermion model also 
found that short-range magnetic order (of the same kind 
that appear in our nematic model) as well as short-range 
orbital order can promote anisotropic transport (59| . 



G. Elastic modulus: ultra-sound measurements 

So far we discussed several manifestations of long-range 
nematic order. An interesting question is how one can 
also probe nematic fluctuations. If these are emergent 
degrees of freedom of the system, then we expect that 
their low-energy excitations will have an important effect 
in the normal state properties. 

As we discussed in Section [III nematic fluctuations are 
four-spin correlation functions, which are rather difficult 
to measure using standard magnetic probes. An alter- 
native is to investigate how they couple to other degrees 
of freedom. For instance, Eq. (JT]) shows that the ne- 
matic order parameter acts as a conjugate field to the or- 
thorhombic distortion. Although in the tetragonal phase 
the mean value of the nematic order parameter is zero, 
nematic fluctuations are still present. As a result, they 
will give rise to a quadratic term in the shear distortion 
e s that effectively reduces the shear modulus C s . Phys- 
ically, nematic fluctuations suppress the energy cost of 
a momentaneous orthorhombic distortion, what is trans- 
lated as a softening of C s . Quantitatively, the renormal- 
ized shear modulus is given by: 



C: 



lim 

<7->0 



d 2 \nZ 
da 2 



(21) 



where Z = J dMi de s e _FmaE_Fcl is the partition func- 
tion, see Eqs. (H} and @. For a harmonic lattice, 
the elastic degrees of freedom can be integrated out 
from the partition function analytically, and we obtain 



Z = J dMi e" 



with F mas given by Eq. 



Taking 



the derivatives with respect to the infinitesimal stress o 
and then making a — > we obtain (27| : 




= c- 1 



\ 2 



\2 

~ s~t Anem 



where Xncm (?) = [{M'( - M| 
fined in Section [TTJ Recall that 



M 2 ) 



(22) 



as de- 
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Figure 9: Typical behavior of the shear modulus in the tetrag- 
onal phase as function of the reduced temperature. The soft- 
ening follows Eq. (|24[) ; here we used k = 0.3. The quantity C s 
is the physical shear modulus, that can be measured in ultra- 
sound experiments, while C 3 is its high-temperature limit. 



IV. NEMATIC ORDER OR ORBITAL ORDER? 

Besides nematic order, it has been proposed that 
ferro-orbital order could drive the tetragonal symmetry- 
breaking observed in the iron pnictides |16h21| . In this 
case, the occupation of the Fe d xz orbital is different 
than the occupation of the d yz orbital. Since the to- 
tal occupation number depends on all states below the 
Fermi level, in some cases it is more convenient to re- 
fer to the splitting between the on-site orbital ener- 
gies, A or b = A xz — A yz ^ 0. ARPES measurements 
have observed A or b < in detwinned electron-doped 
Ba (Fei_ x Co x ) 2 As2 at the same temperature where a 
non-zero resistivity anisotropy sets in. 

By symmetry, the ferro-orbital order parameter enters 
the free energy via: 



rb A 2 
) orb 



M 



(25) 



I q x 2 (g) 
l ~g Lx 2 (q) 



(23) 



with renormalized g = g + \ 2 \/C s (see Eq. [9]). As ex- 
pected, the divergence of the nematic susceptibility leads 
to the vanishing of the shear modulus. 

Most interestingly, Eq. (|22|) shows that one can ex- 
perimentally extract Xnem by just measuring the shear 
modulus C s via standard ultra-sound techniques (27l l60l | . 
To calculate Xnem, one can use inelastic neutron scatter- 
ing data as input, since they provide the parameters of 
X (?) = + q 2 + 7 _1 With this procedure, it is 

possible to establish whether the measured softening of 
C s agrees with what one would expect if the structural 
transition was driven by magnetic fluctuations. This pro- 
cedure was carried on in Ref . (27| , and a very good agree- 
ment was found between the measured C s and the calcu- 
lated Xnem; based on neutron scattering data. 

For a finite-temperature nematic transition, 
I q X 2 (q) ~ ~ (T-T N>0 )-^- d)v . Setting d = 2 

and using the mean- field exponent v = 1/2, we obtain a 
simple expression for the vanishing of the shear modulus: 



C a = CA\ + K 



T-T s 



(24) 



where k is a positive dimensionless constant, of the order 

A 2 

of k ~ c rg-^ , with E mag denoting a characteristic 
magnetic energy scale. Notice that the shear modulus is 
reduced to half of its saturation value at T = T s (l + k). 
The typical behavior of this function is shown in Fig. 
[5] and agrees qualitativel y w ell with the experimentally 
observed softening of C s |27l . [60j | . 



where Xorb = Xorb (q = 0) is the intrinsic ferro-orbital 
susceptibility, with Xorb (?) = (A or b (?) A orb (-?)), and 
A or b is the appropriate coupling constant to the mag- 
netic degrees of freedom. By minimizing the free energy, 
it follows that A or b oc A or b</?- In Ref. [HI, the cou- 
pling constant A or b was calculated using an itinerant ap- 
proach, and was shown to be negative for electron-doped 
compounds, providing a possible explanation for why the 
ARPES measurements observed A or b < when ip > 0. 

Thus, within the nematic scenario, orbital order is in- 
duced by long-range nematic order. However, one could 
argue that in a model with spontaneous orbital order, 
where %~ b ^ T — T or b, orbital order is the one that 
drives the nematic order. Of course, a phenomenological 
approach based only on symmetry considerations cannot 
distinguish between these two different scenarios. One 
then has to rely on model microscopic calculations to 
investigate when spontaneous orbital order and/or spon- 
taneous nematic order appear. 

In what concerns the nematic scenario, as discussed 
in Section [III a spontaneous tendency to long-range ne- 
matic order is found in both strong-coupling (the J\ — 
J2 localized-spin model) and weak-coupling limits (the 
multi-band nested electronic model). On the other hand, 
spontaneous ferro-orbital order is mostly found in par- 
ticular (but not all, see Ref. [lj|) versions of the Kugel- 
Khomskii model, which is obtained in the strong-coupling 
limit. Furthermore, different numerical analyses of the 
five-band Hubbard model did not find a ground state 
that displays orbital order in the absence of long-range 
magnetic order (2^ - |2^ | . Although they do not constitute 
a proof, these results, combined with the fact that the 
magnetism in the iron pnictides is closer to the itiner- 
ant rather than to the localized-spin regime (since these 
systems are metallic, see also optical conductivity mea- 
surements (6l1|). suggest that the nematic instability may 
be the one driving orbital order and structural order. 
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This does not imply that the role played by the orbital 
degrees of freedom is irrelevant: for example, it is clear 
that the reconstruction of the Fermi surface due to or- 
bital order improves the nesting conditions and enhances 
the magnetic transition temperature pjj|62|]. The orbital 
character of the different pockets of the Fermi surface is 
also important for the superconducting state (63| . Nev- 
ertheless, orbital order alone cannot explain the observed 
anisotropic properties of the iron pnictides. For instance, 
calculations of the resistivity anisotropy based on the ex- 
perimentally observed A 0{ h < give the opposite sign of 
the experimental pb — p a |19L 12(1 |22| , indicating that the 
scattering by anisotropic spin fluctuations is fundamental 
to explain the observed resistivity anisotropy. 

Rather than competing tendencies, ferro-orbital and 
nematic order are actually cooperative instabilities, 
which can potentially lead to the enhancement of both 
the structural transition temperature and the anisotropic 
properties of the orthorhombic state. To our knowledge, 
a model that displays both spontaneous nematic and or- 
bital order has not been conceived so far. 

Nevertheless, it is clear that there is a coupling between 
the orbital and magnetic degrees of freedom, since first- 
principle calculations find that the onset of long-range 
SDW order is accompanied by ferro-orbital order, due to 
the different spin polarizations of the d xz and d yz orbitals 
(23 |. If we assume that the system has a tendency to 
orbital order, but that Xnib > at all temperatures, then 
we can integrate out the orbital degrees of freedom from 
the partition function, obtaining the following expression 
for the nematic susceptibility (see Eq. 



Xncm 



K 2 

A. map 



(?) 



— 9 - -"^orb Xorb ^ A el C s 1 (26) 



It is clear from this expression how the different de- 
grees of freedom assist nematic fluctuations and nematic 
order. The first term is just the usual Gaussian fluctua- 
tions of the magnetic order parameter, going to zero only 
at the (bare) magnetic transition temperature. The sec- 
ond term g > also comes from the magnetic sector (see 
Eq. [IJ and is responsible for the degeneracy of the two 
magnetic stripe states, i.e. it is the "intrinsic" magnetic 
instability towards a nematic state. The third term is 
associated with the orbital degrees of freedom and shows 
that a tendency towards ferro-orbital order enhances the 
nematic susceptibility. Finally, the last term, which orig- 
inates from the magneto-elastic coupling, implies that 
softer lattices have a larger effective nematic coupling. 



V. INTERPLAY BETWEEN NEMATICITY 
AND SUPERCONDUCTIVITY 

In the previous sections, we explored the effects of the 
nematic degrees of freedom on the normal state prop- 
erties. Here, we investigate how superconductivity and 
nematic order/fluctuations affect each other. 



A. Competition with superconducting order 

It has been established both experimentally and theo- 
retically that SC and SDW are competing orders in the 
iron pnictides. Neutron diffraction measurements ob- 
serve a strong suppression of the SDW order parame- 
ter below T c - in some materials, there is even a reen- 
trance of the non-magnetic phase at low temperatures 
Q- Theoretically, the fact that the same electrons are 
responsible for both the static staggered magnetization 
and the superfluid density explains such a strong com- 
petition [6J, [6a|. In our phenomenological approach, we 
can capture this competition by adding to the free energy 
expansion a bi-quadratic SDW-SC term: 



MS) A 



F [Mi, A] = F mag [Mi] + F sc [A] + A sc (M 2 

(27) 

In the mean-field level, a non-zero SC order parameter 
renormalizes the magnetic susceptibility, suppressing the 
SDW transition temperature Tn 0' 



A S ,A^ 



(28) 



Within this phenomenological approach, it is straight- 
forward to understand qualitatively the effects of super- 
conducting order on the nematic order parameter. Below 
T c , the magnetic susceptibility is suppressed according to 
Eq. (|28p . Since long-range nematic order is induced by 
magnetic fluctuations, it will also be suppressed by SC. 
This implies that both states are also competing orders, 
despite the fact that there is no direct coupling between 
the nematic and SC order parameters. 

To investigate the interplay between nematicity and 
SC quantitatively, we consider the case where the system 
has no long-range magnetic order once SC sets in. First 
we analyze the behavior of the nematic order parameter 
ip, i.e. we assume T c < T s . For simplicity, we solve 
the self-consistent equations ([2]) in d = 2, since in this 
case there is no finite-temperature SDW transition. A 
straightforward calculation leads to an implicit equation 
for dp = 4:nip/gT N (see Ref. (U): 



u.( <p 

w coth <p H — In 

g \ sinh dp 



r + A SC A 2 



(29) 



where fo = a (T — T/v) and A sc = Air\ sc /gTN- For a 
fixed value of the ratio u/g, it is straightforward to solve 
the self-consistent equation (|29|) for (p. In mean-field, for 
T > T c , we have A 2 = 0, whereas for T < T c , it holds 
that A 2 = Aq (T c — T) . The structural transition tem- 
perature takes place at T s = Xjv + a -1 . Using Eq. (|2"9"| . 
we can evaluate dp/dT for all temperatures, obtaining: 



dip 
dT 



T C <T<T S 
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Figure 10: Nematic order parameter ip (in units of its zero- 
temperature value ipo) as a function of temperature (in units 
of the structural transition temperature T s ). The dashed line 
corresponds to the case with no SC, while the solid line corre- 
sponds to the solution in the presence of a SC state that sets 
in at T c = T a /2. The SDW-SC coupling constant was set to 

= o.5 (a), 1.5 (b) and 2.3 (c). 



dip 
dT 



(30) 



T<T C 



with c > a temperature-dependent positive number. 
Eqs. (f30|) show that the nematic order parameter is sup- 
pressed below T c . If the competition between SDW and 
SC is strong enough, such that A sc Aq > a, then ^ 
changes sign below T c , implying that the orthorhombic 
distortion actually decreases inside the SC state. This 
opens up the possibility of completely killing the Ort 
phase at low temperatures, promoting a reentrance of 
the Tet phase. In Fig. [TU1 we illustrate these behaviors 
by plotting the solutions of Eq. (|29|) as a function of 
temperature for different values of the coupling constant 
A sc . The competition between the two order parameters 
has also been investigated by Ref. (66|, where a different 
approach was employed. 

Remarkably, x-ray diffraction measurements have 
found a very strong suppression of the orthorhombic dis- 
tortion inside the SC phase, and even a reentrance of 
the tetragonal phase for Ba (Fci_ x Co x ) 2 As2 near opti- 
mal doping Q. This behavior was observed in the ab- 
sence of long-range magnetic order, when only the SC 
and nematic order parameters are non-zero. These mea- 
surements provide yet another strong evidence for the 
electronic character of the structural transition in the 
iron pnictides, and find a natural explanation within the 
nematic scenario. No additional coupling between the 
nematic and SC order parameters is necessary, as the 
magnetic origin of the nematic phase, combined with the 
competition between SDW and SC, give rise to an indi- 




Figure 11: Shear modulus C 3 (in units of its high-temperature 
value C s ) as a function of temperature (in units of T c ). The 
dashed line refers to the case without SC and with a putative 
SDW quantum critical point. The solid line corresponds to 
the case where SC sets in at T c , and reveals the hardening of 
the shear modulus inside the Tet-SC phase. 



rect competition between nematicity and SC. 

This competition is also manifested in the shear modu- 
lus C s . An interesting case to analyze is when the system 
does not have long-range nematic or magnetic order, but 
only SC at a finite T c (as it is the case in optimally doped 
samples). For concreteness, we consider the presence of a 
SDW quantum critical point, such that rocT. To calcu- 
late the shear modulus, we need to evaluate J Xmag (l) m 
Eq. (|23|l . Above T c , we have the usual overdamped spin 



excitations, x 



mag 



(?) 



V +7 



simplicity d = 2 and the T 



u> n | . Considering for 
limit, we obtain: 



Xnom 



47T Z V T7 



(31) 



where A is the frequency cutoff. Below T c , there are 
changes not only in the static part of the susceptibility, 
r — > r = r + A SC A 2 , but also in the dynamics of Xmag (<?)■ 
The presence of a SC gap makes the spin dynamics ballis- 
tic for energies w„ < A, while for w„ ^> A one recovers 
the overdamped dynamics (67| . To capture these two 
different behaviors, we write the phenomenological form 
Xmag (?) = r + q 2 + f (u> n ) with / (ui n ) = u?J (A7) for 
u> n < A and / (uj n ) = 0J n /^ for tu n > A. A straightfor- 
ward calculation then gives: 



Anem 



47T 2 
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with f = r + A SC A 2 . Substituting this expression in Eq. 
(|22p . we obtain the renormalized shear modulus. In Fig. 
ITT| we plot C s using parameters such that 7 3> A, in 
accordance to what neutron scattering experiments ob- 
serve. Notice that the shear modulus is enhanced below 



T c , i.e. the lattice is hardened in the SC phase, due to the 
suppression of nematic fluctuations. We emphasize that 
the origin of this behavior is in the magnetic nature of the 
nematic degrees of freedom, allied to the competition be- 
tween SDW and SC. The hardening of the shear modulus 
below T c has been observed by ultra-sound measurements 
in optimally doped Ba (Fci_ x Co x ) 2 As2 [27l. |60|. provid- 
ing another evidence in favor of the nematic model. 

B. Nematic fluctuations and superconductivity 

Given the evidence in favor of nematic degrees of free- 
dom in the iron pnictides, it is natural to ask whether 
nematic fluctuations are able to enhance the value of T c . 
Here we give a qualitative argument that this may be the 
case in the iron pnictides, although more rigorous calcu- 
lations are necessary to confirm it. We start with a sim- 
plified model containing one hole pocket and one electron 
pocket, both with the same density of states Np. The 
linearized BCS equations for this two-band problem is 
given by the matrix equation [68| : 

it)"** (33) 

where W is the typical energy scale of the pairing inter- 
action, Ui and U2 are intra-band interactions, and V is 
the inter-band interaction. The value of T c is given by 
the largest eigenvalue A+ of the matrix: 

T c = We- 1/iNpX+) (34) 

If Ui and U2 are repulsive interactions, i.e. Ui > 0, 
then SC will only appear if \V\ > \jU\Ui- The largest 
eigenvalue is given by: 

^ = -(^)+\/ffl 2 + - (35, 

and the eigenvectors give Ai and A2 with the same sign 
if V < (attractive) or Ai and A2 with opposite signs 
if V > (repulsive). In the first case, the SC state is 
called s ++ state, while in the second case it is called s H 
state. In either case, the interband interaction must be 
large enough to overcome the intraband repulsion. In 
the s H state this becomes possible due to the proxim- 
ity to a SDW instability, which enhances the interband 
repulsion V 0|. The s ++ state is the leading instability 
if interband orbital fluctuations (i.e. anti-ferro orbital 
fluctuations) are strong enough to enhance the electron- 
phonon attractive interaction (69j . Several experimental 
results indicate that the state realized in the iron pnic- 
tides is the s H . Evidence includes the existence of a 
resonance mode in the magnetic excitation spectrum be- 
low T c (70| ; the pattern of quasi-particle interference in 
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the presence of a magnetic field (71) : half flux-quantum 
transitions in composite loops of Nb-iron pnictides [72| ; 
and the microscopic coexistence between SC and SDW 
1- 

How nematic fluctuations affect this model? The key 
point is that, unlike the SDW fluctuations, nematic fluc- 
tuations are peaked at q = 0. Therefore, their main 
effect should be on the intraband component of the in- 
teraction. Moreover, unlike magnetic fluctuations, ne- 
matic fluctuations couple to the charge of the electrons. 
As such, their coupling the electrons is analogous to the 
coupling between the latter and acoustic phonons, result- 
ing in an attractive intraband interaction U ncm < 0, i.e. 
U\ — > U\ — |f/ nem |. Then, according to Eq. (|35|) . the 
largest eigenvalue A + would increase and, consequently, 
T c would be enhanced. Notice that since the pairing in- 
teraction due to nematic fluctuations has an intra-band 
character, it affects T c in the same way regardless of 
whether the ground state is s ++ or s H . 

Of course, further calculations are necessary to confirm 
this qualitative picture. One argument that supports a 
rather weak additional role of nematic-pairing is the ab- 
sence of significant variations in the gap amplitude on the 
hole pocket (although recent measurements in LiFeAs 
found angular variations in the hole-pockets gaps, see 
[ill]). Notice that the nematic fluctuations would always 
couple to the form factor (fc^ — k^j , which should give 
rise to smaller pairing gaps on the Brillouin zone diago- 
nals. Also, key points neglected in this analysis are the 
possibility of different frequency cutoffs W among the dif- 
ferent interactions and the fact that the nematic coupling 
should be treated as a quantum critical pairing where all 
energy scales of the problem contribute to the pairing in- 
teraction, not only the ones that are below the typical 
excitation energy of the nematic mode (74|. Yet, since 
nematic fluctuations are generated by magnetic fluctua- 
tions, it is certainly an interesting open problem to in- 
vestigate whether this novel collective mode can play a 
favorable role in Cooper pairing. 



VI. CONCLUDING REMARKS 

In this paper we discussed distinctive manifestations 
of the nematic degrees of freedom in several properties of 
the iron pnictides. Due to the magnetic origin of the ne- 
matic phase, magnetic properties such as the spin-lattice 
relaxation rate, the spin-spin correlation function, and 
the uniform magnetic susceptibility display characteristic 
features when the nematic order parameter becomes non- 
zero at T s . Furthermore, the enhancement of magnetic 
fluctuations at T s also has consequences to the electronic 
spectrum, opening a pseudogap in the electronic spectral 
function. Scattering of electrons by these anisotropic spin 
fluctuations in the nematic paramagnetic phase leads to 
a resistivity anisotropy whose sign is different in electron- 
doped and hole-doped compounds. 

We also showed that the magneto-elastic coupling 
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makes the orthorhombic distortion proportional to the 
nematic order parameter, and the shear modulus propor- 
tional to the inverse nematic susceptibility. In the normal 
state, the shear modulus is softened by low-energy ne- 
matic fluctuations, but in the SC state it becomes harder 
due to the indirect competition between nematicity and 
SC. This indirect competition, rooted on the direct com- 
petition between SDW and SC for the same electronic 
states, is also manifested in the suppression of the or- 
thorhombic distortion below T c . The role of orbital or- 
der to the onset of long-range nematic order was also 
discussed. We showed that even when orbital order is 
not a spontaneous instability of the system, ferro-orbital 
fluctuations enhance the nematic susceptibility, whose di- 
vergence then gives rise to a finite orbital polarization. 

It is important to emphasize that, besides orbital order 
and nematic order, other mechanisms have been proposed 
for the electronic tetragonal symmetry-breaking in the 
pnictides. In Ref. (75|, the authors argue that an orbital- 
current density wave may be spontaneously generated 
due to the nesting features of the Fermi surface. As a 
secondary effect of this orbital-current state, a charge- 
density wave appears, triggering the structural transi- 



tion. The authors of Refs. [76|, [77|, on the other hand, 
suggest proximity to an orbital-selective Mott transition 
as the origin of the orthorhombic state. Kontani et al. 
propose an electronic nematic state rooted on quadrupo- 
lar orbital interactions allied to impurity effects [78[. In 
Ref. (79^ . Daghofer et al. performed a numerical study 
of a three-band Hubbard model and found nematic fea- 
tures in the electronic spectral function, although the 
tetragonal symmetry of the model was explicitly bro- 
ken by an anisotropic exchange constant. In Ref. (80j |. 
Paul proposes a mechanism for the structural transition 
based on the magneto-elastic coupling, which is very sim- 
ilar to the emergent nematic model discussed here. We 
note also that some works treat the coupling between the 
SDW and structural transitions phenomenologically [8l| . 
Such an approach has the clear advantage of providing 
model-independent predictions, but it does not address 
why these two transitions follow each other in the phase 
diagram. 



Several of the manifestations of the nematic phase dis- 
cussed here have been detected experimentally - such as 
the suppression of the orthorhombic distortion below T c 
(x-ray diffraction) the enhancement of 1/XiT below 
T s (NMR) [48[ ; the anisotropy in the uniform suscep- 
tibility (torque magnetometry) [I^; the softening of the 
shear modulus in the normal- Tet phase and its hardening 
in the SC-Tet phase (ultrasound measurements) (27l . [60j ; 
and the change in the sign of the resistivity anisotropy 
upon hole-doping (transport) [l5T| . Yet, there are two im- 
portant predictions of the nematic model that remain to 
be seen: the inequivalence between magnetic fluctuations 
around the ordering vectors Qi = (n, 0) and Q 2 = (0, it) 
below T s (inelastic neutron scattering) and the possible 
opening of a pseudogap in the electronic spectrum below 
T s (ARPES and STM). Other manifestations of nematic 
order not presented here have also been discussed in the 
literature, such as the effects of nematicity on the SC 
vortex structure [12, 8.2], on the elastic domain walls [83^ . 
and on the scattering by impurities Q. 

Overall, the nematic model offers a simple framework 
to understand the interplay between the several degrees 
of freedom present in the phase diagrams of the iron 
pnictides, shedding light on the primary role played by 
magnetism. An interesting topic that deserves further 
attention is the importance of nematic fluctuations to 
high-temperature SC. Here, we gave a qualitative argu- 
ment of how T c can be enhanced by nematic fluctua- 
tions. A detailed investigation of this result can poten- 
tially have implications also for other systems where ne- 
maticity has been proposed, such as some cuprates [HJ 
and some heavy fermion systems (85| . 
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